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ABSTRACT 

Recently, the temperature T and luminosity Lx of the hot gas halos of early type galaxies have been 
derived with unprecedented accuracy from Chandra data, for a sample of 30 galaxies covering a wider 
range of galactic luminosity (and central velocity dispersion <Tc) than before. This work investigates 
the origin of the observed temperatures, by examining the relationship between them and the galaxy 
structure, the gas heating due to Type la supernovae (SNIa's) and the gravitational potential, and 
the dynamical status of the gas flow. In galaxies with (Tc^200 km s~^, the T's are close to a fiducial 
average temperature for the gas when in outflow; at 200< crc(km s~^) < 250, the T's are generally 
lower than this, and unrelated with ac, which requires a more complex gas flow status; at larger ctc, 
the T's may increase as a^, as expected for infall heating, though heating from SNIa's, independent of 
(Tc, should be dominant. All observed T's are larger than the virial temperature, by up to ~ 0.5 keV. 
This additional heating can be provided in the X-ray brightest galaxies by SNIa's and infall heating, 
with a SNIa's energy input even lower than in standard assumptions; in the X-ray fainter ones it 
can be provided by SNIa's, whose energy input would be required close to the full standard value at 
the largest CTc- This same energy input, though, would produce temperatures larger than observed 
at low ac, if entirely thermalized. The values of the observed T's increase from outflows to inflows; 
the gas is relatively hotter in outflows, though, if the T's are rescaled by the virial temperature. For 
200 < (Tc(kin s^^) < 250, lower Lx values tend to correspond to lower T's, which deserves further 
investigation. 

Subject headings: galaxies: elliptical and lenticular, CD — galaxies: fundamental parameters — 
galaxies: ISM — galaxies: kinematics and dynamics — X-rays: galaxies ~ X-rays: 
ISM 



^ . 1. INTRODUCTION 

OO I The advent of the Chandra X-ray observatory, with 
Q\ . its unprecedented sub-arcsecond resolution, allowed to 
OO study better than ever before the main contributors to 
CN the total X -ray emission of early-type galaxies (hereafter 
\^ ETGs): the low-mass X-ray binaries (LMXBs; Fabbiano 
2006), a population of weak sources as late type stel- 
7-H lar coronae, cataclismic variables, and coronally active 
I , binaries (Pellegrini & Fabbiano 1994, Revnivtsev et al. 
;v ' 2008), the nuclear emission due to a supermassive black 
hole (MBH; e.g., Gallo et al. 2010, Pellegrini 2010), and 
K> , a hot interstellar medium (ISM) with a temperature of a 
5_] ■ few million degrees. After careful subtraction of the stel- 
es ■ lar (resolved and unresolved) and nuclear emissions, the 
properties of the hot ISM could be characterized with 
unprecedented accuracy. Recently, this has been done 
for a sample of 30 normal (non-cD) ETGs observed with 
Chandra to a depth ensuring the detection of bright 
LMXBs (Boroson et al. 2011, hereafter BKF). This is 
the first X-ray sample of ETGs covering a wide range 
of galactic luminosity, central velocity dispersion ctc, and 
hot gas emission Lx, and with the X-ray properties of 
the hot gas (e.g., luminosity Lx and average temperature 
T) derived in a homogeneous way, using a complete and 
accurate procedure to subtract all kinds of non-gaseous 
emission (nucleus, detected and undetected LMXBs, and 
unresolved weak stellar sources). This approach resulted 
in a larger fraction of hot gas-poor galaxies than in pre- 
vious samples, with Lx extending down to much lower 
values than before (~ 10'^^ erg s~^), and showing a varia- 



tion of up to ~ 3 orders of magnitude at the same galactic 
luminosity (see also David et al. 2006, Diehl & Statler 
2007, Memola et al. 2009). Such a wide variation, even 
larger than previously found, had been linked to the ori- 
gin and evolution of the hot ISM, and had provided ev- 
idence for the effectiveness of an internal heating mech- 
anism (as from type la supernovae, hereafter SNIa's) to 
regulate the gas evolution and produce its very differ- 
ent content in ETGs at the present epoch (Loewenstein 
& Mathews 1987, David et al. 1990, Ciotti et al. 1991); 
the action of external agents (gas stripping, confinement, 
accretion) to reduce or enhance the gas content was also 
invoked (e.g.. White & Sarazin 1991, Brown & Bregman 
2000, Sun et al. 2007). 

With this new characterization of the hot gas, BKF re- 
visited the relationships between fundamental properties 
of the hot gas and of the host galaxy, as the Lx — T, 
Lx — (Jc and T — ac relations, where CTc is a representa- 
tive measure of the depth of the galactic potential well 
(Eskridge et al. 1995, O'Sullivan et al. 2001, 2003). Lx 
correlates positively with T and ac, though with a wide 
variation at fixed <7c and T. Interestingly, the best fit re- 
lation Lx oc T**'^, close to what already known for X-ray 
luminous ETGs (O'Sullivan et al. 2003), is still moder- 
ately strong among ETGs with low T and Lx', also, in the 
Lx - ac relation, ETGs with kT > 0.4 keV are the X-ray 
brightest (with one exception), while those with kT < 0.3 
keV are the X-ray faintest. The least gas rich ETGs are 
then the coolest ones, which seemed contrary to expecta- 
tion, if low Lx ETGs loose their ISM in an outfiow (e.g.. 
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David et al. 1990, Ciotti et al. 1991), and the hotter the 
gas, the stronger is the outflow (BKF). On average T 
increases with ac, and most ETGs lie above a rough esti- 
mate of the gas virial temperature (To- = ^iripa'^/k), sug- 
gesting the presence of additional heating. ETGs with a 
moderate to high gas content {Lx > 5 x 10'^^ erg s~^) 
follow a trend roughly parallel to that of T„; instead, 
ETGs with little hot gas {Lx < 5 x 10^^ erg s"^) have a 
similar temperature for CTc ranging from 160 to 250 km 
s~^. This lack of correlation was attributed to a different 
dynamical state of the hot ISM in gas-poor with respect 
to gas-rich ETGs, though a full explanation of this aspect 
remained to be found (BKF). 

This work takes advantage of the new accurate mea- 
surements of the hot gas properties, and of the fundamen- 
tal relations Lx — ctc and T — derived down to galaxy 
masses and X-ray luminosities smaller than ever before 
(BKF), to investigate the relationship between T, the 
galaxy structure, the internal gas heating mechanisms 
(SNIa's, and those linked to the gravitational potential), 
and the dynamical status of the gas flow. To this pur- 
pose, a few characteristic temperatures are introduced, 
depending on the nature of the gas heating sources and 
the galaxy structure, and relevant for the various gas flow 
phases; these characteristic temperatures are then com- 
pared with the observed T values. In doing so, galaxy 
mass models are built according to the most recent un- 
derstanding of the ETGs' structure, such as their stellar 
mass profile and their dark matter content and distri- 
bution, as indicated by detailed modeling of optical ob- 
servations and by the main scaling laws (e.g., Cappellari 
et al. 2006, Weijmans et al. 2009, Auger et al. 2010, 
Napolitano et al. 2010, Shen & Gebhardt 2010). The 
aims are to address the following questions: can the gas 
heating sources above account for the observed T's? how 
are the various input energy sources for the gas used in 
the different flow phases? is there any relation between 
T and the flow phase? 

We present in Sect. [2] the sources of mass and heating 
for the hot ISM, in Sect. [3] the conditions for the gas 
to escape from the galaxy, in Sect. |4] the galaxy mass 
models, in Sect. [5] the comparison between observed and 
predicted temperatures, in Sect. IH] the relation between 
gas temperature and flow status, and in Sect. [7] the con- 
clusions. 

2. SOURCES OF MASS AND HEATING FOR THE HOT GAS 
2.1. Gas Mass 

In ETGs the hot gas comes from stellar mass losses 
produced by evolved stars, mainly during the red giant, 
asymptotic giant branch, and planetary nebula phases, 
and by SNIa's, that are the only ones observed in an old 
stellar population (e.g., Cappellaro et al. 1999). The 
first, more quiescent, type of losses originates ejecta that 
initially have the velocity of the parent star, then individ- 
ually interact with the mass lost from other stars or with 
the hot ISM, and mix with it (Mathews 1990, Parriott 
& Bregman 2008). 

For a galaxy of total stellar mass M*, the evolution of 
the stellar mass loss rate M*(t) can be calculated using 
single burst stellar population synthesis models (Maras- 
ton 2005), for a Salpeter and for a Kroupa Initial Mass 
Function (IMF), assuming for example solar abundance. 



So doing, at an age of 12 Gyrs, a rate is recovered of M,= 
B Lb{Lb^q) Moyr~^, where Lb is the galactic 

B-band luminosity at an age of 12 Gyr, and B=1.8 or 
B=1.9 for the Salpeter or Kroupa IMF (see also Pellegrini 
2011). This value is in reasonable agreement with the av- 
erage derived for nine local ETGs from ISO data (Athey 
et al. 2002) of = 7.8 x lO^^^ Lb{Lb,q) Moyr-\ an 
estimate based on individual observed values that vary 
by a factor of ^ 10, though, which was attributed to 
different ages and metallicities. 

The total mass loss rate of a stellar population M is 
given by the sum M = Af, + Msn, where Msn is the rate 
of mass loss via SNIa events for the whole galaxy. Msn 
is given by Msn = MsnRsn, where Msn = 1.4M0 is 
the ejected mass by one event, and i?sN is the explosion 
rate. i?sN has been determined for local ETGs to be 
i?SN - 0.16(i?o/70)2 X lO-^^LB{LB,Q)yT^\ where Ho 
is the Hubble constant in units of km s~^ Mpc^^ (Cap- 
pellaro et al. 1999). More recent measurements of the 
observed rates of supernovae in the local Universe (Li et 
al. 2010) give a SNIa's rate in ETGs consistent with that 
of Cappellaro et al. (1999). For Ha ^ 70 km s^^ Mpc"\ 
one obtains Msn =2.2x 10""Ls(Lb,o) Mq yr"\ that 
is ^ 80 times smaller than M* derived above for an age of 
12 Gyr; therefore, the main source of mass for the hot gas 
is provided by M*. A reasonable assumption is that the 
gas is shed by stars with a radial dependence that follows 
that of the stellar distribution, so that the density profile 
of the injected gas is Pgasir) oc p*{r), where p*{r) is the 
stellar density profile. This assumption is adopted here- 
after, and the characteristic temperatures presented be- 
low apply to a gas distribution following pgas (r) oc p^, (r) 
(but see also Sect. [5]). 

2.2. Heating from stellar motions and supernovae 

The material lost by stars is ejected at a velocity of 
few tens of km s~^ and at a temperature of <10^ K 
(Parriott & Bregman 2008), and is subsequently heated 
to high. X-ray emitting temperatures by the thermaliza- 
tion of the stellar velocity dispersion, as it collides with 
the mass lost from other stars, or with the ambient hot 
gas, and is shocked. Another source of heating for the 
stellar mass losses is provided by the thermalization of 
the kinetic energy of SNIa's events. The internal en- 
ergy given by these heating processes to the unit mass of 
injected gas is 3kTinj /2pmp (with k the Boltzmann con- 
stant, nip the proton mass, prup the mean particle mass, 
with p — 0.62 for solar abundance); Tinj is determined 
by the heating due to thermalization of the motions of 
the gas- losing stars {Tstar), and of the velocity of the 
SNIa's ejecta (Tsn), and is written as (e.g., Gisler 1976, 
White & Chevalier 1983): 



T 

-' 1.7 



M,T, + MsnT, 



SN 



M 



(1) 



Here is the equivalent temperature of the stellar mo- 
tions (see below), and Tej — 2pmpEsN / {3kAIsN) is the 
equivalent temperature of the kinetic energy i?sN of the 
SNIa's ejecta, with Esn = 10^^ erg for one event (e.g., 
Larson 1974). T^j can be calculated assuming that a 
factor / of is turned into heat; / < 1, since radia- 
tive energy losses from expanding supernova remnants 
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may be important, and values down to / = 0.1 have 
been adopted (Larson 1974, Chevalier 1974); a value of 
/ = 0.85 could be not too far off for the hot diluted 
ISM of ETGs (e.g., Tang & Wang 2005). In this way, 
Tej = (//0.85)1.5x lO^K. From approximating M ~ M*, 
and using the estimates of Sect. 12.11 for Msn and 
(for a Kroupa IMF) at the present epoch, one obtains 
TsN ^ 1.7(//0.85) X 10^ K. 

The injection temperature Tinj is then the sum of two 
parts: one (Tsn) is independent of the position within 
the galaxy where the gas is injected (e.g., independent 
of radius in spherical symmetry), and is also constant 
from galaxy to galaxy (for fixed IMF and age of the stel- 
lar population, and SNIa's rate); for each ETG it can, 
though, evolve with time, if Mg^ and M* evolve dif- 
ferently with time (Ciotti et al. 1991). The other part 
{Tstar) is instead basically independent of time, but has 
a radial dependence, and changes with the galaxy struc- 
ture, i.e., with the total mass and its distribution. An 
average T* is obtained calculating the gas mass- weighted 
temperature gained by the thermalization of the stellar 
random motions, < >: 

<T*>-1^ Jinr^Mr)a^ir)dr, (2) 

where cr(r) is the one-dimensional velocity dispersion of 
the stars. The integral term in Eq. [2] is the same that 
gives the kinetic energy associated with the stellar ran- 
dom motions [Ekin = 1-5 / inr-^ p^,{r)a^ (r) dr], and that 
enters the virial theorem for the stellar component; the 
mass- weighted temperature in Eq. [2] is then often called 
"gas virial temperature" . For a galaxy mass model made 
of stars and dark matter, characterized hy TZ — Mh/M^,, 
where Mh is the total dark mass, and /3 — rh/r,,, with 
r/i and the scale radii of the two mass distributions, 
< T* > can be expressed using the central velocity dis- 
persion CTc as < T, fj,nipa'^V,{TZ, l3)/k (e.g., Ciotti 
& Pellegrini 1992). The function ft increases mildly for 
larger TZ and for lower /3, that is for a larger amount of 
gravitating mass or a higher mass concentration, but al- 
ways < 1, since a(r) has in general a negative radial 
gradient (e.g.. Sect. |3] and Fig. [1] below). < T, > is then 
proportional to , and a simplified version of the virial 
temperature in Eq.[2]that is often used is = /^TOpCr^/fc; 
To- of course overestimates the true < > . 

The mass-averaged injection temperature is finally 
given by 

< T,nj >=< n > +1.7(//0.85) X 10^ K, (3) 

where in general the second term dominates, as is shown 
in Sect. [5] below. 

2.3. Heating during Infall 

In case of mass losses flowing to the galactic center, 
the gas can be heated due to infall in the galactic po- 
tential and adiabatic compression; this process is some- 
times referred to as "gravitational heating" . The average 
change in gravitational energy per unit gas mass inflow- 
ing through the galactic potential down to the galactic 
center is 

= TT / 47rrV.(r)[0(r) - (4) 



for galaxy mass distributions with a finite value of 0(0) 
(see also Ciotti et al. 1991). One can define a temper- 
ature equivalent to the energy in Eq. |4] as < T^,.^^ >= 

2^mp E^^^^/3k. As < T, >, also < T^^^ > is cx tr^, and 
increases for larger TZ and smaller /3, which, for inflowing 
gas, can be understood as a larger gas heating by com- 
pression during infall for a larger dark matter amount or 
its higher concentration. 

Not all of i?^.av can be available for heating, though. 
If the inflow keeps quasi-hydrostatic, then, by the virial 
theorem, the energy radiated away is roughly one-half 
of the change in the gravitational potential energy, and 
that available for the heating of the gas is the remaining 
half] (i.e., ^ O.dE^j.^^). Actually, the energy available for 
heating will be much less than this. Inflows are caused by 
the radiative losses produced by the accumulation of the 
stellar mass return, that makes the cooling time lower 
than the galactic age; in the central regions, within a 
radius of ^ 1 kpc, the cooling time can be as short as 
<10® yr, even shorter than the infall time (e.g., Sarazin 
& White 1988, Pellegrini 2011). In these conditions, the 
gas departs from a slow inflow, becomes very dense and 
supersonic close to the center, and cools rapidly down to 
low temperatures, so that > O.bE^^.^^ is radiated away 
or goes into kinetic energy of condensations (Sarazin & 
Ashe 1989). Furthermore, there is the possibility that 
not all the gas reaches hot the galactic central region, if 
thermal instabilities develop and produce drop-outs from 
the flow; if gas cools and condenses out of the flow at 
large radii, then -E^av can be much lower than in the 
deflnition above, and heating due to infall in the grav- 
itational potential is "lost" (Sarazin & Ashe 1989). In 
conclusion, without a precise knowledge of how to com- 
pute -B^-av (which depends on the radius at which the 
injected gas drops below X-ray emitting temperatures), 
and about what fraction of iS^av is radiated or goes into 
kinetic energy of the condensations, < T+a^ > remains 
a reference value; a more direct use can instead be made 
of the analogous temperature for escape < T~^^ > intro- 
duced in Sect. [3] below. 

2.4. Heating from a central MBH 

Due to the presence of a central MBH in ETGs, an- 
other potential source of heating for the gas could be 
provided by nuclear accretion. This subject has been 
studied intensely recently, through both observations and 
modeling, and it appears that the energy provided by ac- 
cretion is of the order of that needed to offset cyclically 
the cooling of the inflowing gas in the central regions of 
gas-rich ETGs (e.g., Birzan et al. 2004, MiUion et al. 
2010). Therefore, the central MBH is believed to be a 
heat source that balances the radiative losses of the gas, 
acting mostly in the central cooling region. In gas poor 
ETGs the nuclear accretion energy is far lower due to 
the very small mass accretion rate, if present (Pellegrini 
et al. 2007), and the absorption of the energy output 
from accretion is likely not efficient. Given the role of 
the MBH outlined above, possible energy input from the 

^ The energy lost in radiation and that converted into heat are 
actually each equal to 0.5-E^av for a self-gravitating gas; for gas 
in an external potential, the result should remain roughly valid 
(Binney & Tremaine 1987). 
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MBH will not be considered as a source of global heating 
for the gas. 

3. CONDITIONS FOR ESCAPE 

Another characteristic temperature for comparison 
with observed T values is the temperature with which 
the gas can escape from the galaxy. Assuming that 
the flow is stationary and adiabatic, the Bernoulli con- 
stant on each streamline along which the gas flows out 
of the galaxy must be positive. The Bernoulli equa- 
tion with the minimum energy for escape is written as 

H{r) + v^{r)l2 + 4>{r) = 0, where i/ ^ = ^ is 

the enthalpy per unit gas mass, 7 is the ratio of specific 
heats, Cs is the sound velocity, and v is the flow veloc- 
ity. Integrating over the galaxy volume and gas-mass 
averaging, this condition becomes: 

/ 4TTr'^p^{r)H{r)dr+- Aur'^ p^{r)v'^{r)dr = M^E'.^ 
Jo ^ Jo 

(5) 

where the escape energy 

1 f°° 

E;... = -JfJ^ ^7rr'p*{r)Hr)dr (6) 

is the average energy required to remove from the galaxy 
the unit gas mass. E~^g^ gives a minimum energy re- 
quirement, since energy losses due to cooling may be 
present; but these are not important for outfiows that 
typically have a low density (e.g.. Sect. 15.31 below). The 
escape temperature equivalent to £'~av ^ ^grav > = 
2iimp E~j.^^/3k. The condition for the minimum energy 
for escape can then be translated into a condition for the 
injection temperature of the gas < Ti„j > to be larger 
than < T~^^ >. In the simple case that 0(r) is due 
only to one (stellar) mass component, from < T, > = 
2/impiJfci„/3fcM* and the virial theorem (Sect. [5]), one 
derives that < T^^.^^ 4 < T, >. The estimates of 
E~^^ in Eq. [6] and then of < T^^^^ > will be calculated 
below for a general mass model (e.g., made by the super- 
position of stars and dark matter, with a different radial 
distribution). In previous works, the sufficient condition 
for the existence of a galactic wind was that the injection 
temperature exceeded an "escape temperature" , defined 
as 2To- (White & Chevalier 1983), or "twice the equiv- 
alent dark halo temperature" (Loewenstein & Mathews 
1987), coupled with the request for the radiative cooling 
time in the central part of the galaxy to be longer than 
the time required to flow out of this region. These con- 
ditions are similar to imposing that < Tinj > exceeds 
< Tg^av > derived above. 

In principle the gas can escape with different combina- 
tions of V and T, and the observed T should be closeO to 
that entering H in Eq. [S] There are two extreme cases for 
the value of the flow velocity v with respect to Cg (i.e., to 
the temperature). One is when the material is brought 
to infinity keeping a subsonic velocity, then the minimum 

^ For example, we recall two approximations made here with 
respect to the case of real ETGs: the total gas profile may be 
different from that of the stars, and the flow has a time-continuous 
distributed mass and energy input. The first of these points will 
be further discussed in Sect. |5] 



energy requirement becomes iJ « — (/>; neglecting the ki- 
netic term in Eq. [5l one then obtains a characteristic 
gas-mass averaged subsonic escape temperature: 

< If" r ^^"-'p*^'^^^'^^' = I < ^^-v > . 

/7) 

In a general case, H ^ —cf) gives for T a larger require- 
ment than obtained when v is not neglected; also the 
partition of the gas energy between enthalpy and kinetic 
energy can vary with radius within a galaxy. All this 
means that < T^^c > represents a fiducial upper limit 
to the observed temperatures of outflowing gas: if the 
kinetic energy of the flow is important, then the actual 
gas temperature will be lower (the stronger the outflow, 
with respect to Cg, the cooler the gas). At the opposite 
extreme case where the "temperature" contribution to 
the gas energy is minor and that of the velocity is dom- 
inant, the Bernoulli equation reduces to v\^j2 -I- (/> = 0; 
this gives the usual escape velocity of a unit mass from 
a potential well: Wesc('') = •\/2|(/)(r)|. 

Finally, < T-^._^^ > and < T^i^'J' > have the same de- 
pendence as < > on (7^,7?., /?. For the represen- 
tative galaxy mass models used here (Sect. SI Fig. [2]), 
^grav — (1-^ — 2.3)ii^~av; since <0.5i?+av can be con- 
verted into heat, then the corresponding temperature 
gained by infaU wfll be Tinfi< < T'g7av >■ 

4. GALAXY MASS MODELS 

In this work the exact values of < >, < T^^^ > and 

< Tgi^av > calculated as a function of CTc, for a series 
of representative 3-component galaxy mass models, made 
by the superposition of a stellar distribution and a dark 
matter halo, to which a central MBH is added. The stel- 
lar density profile is given by the deprojection of a Sersic 
law with index n = 4 or 5, as appropriate for ETGs of 
the luminosities considered in this work (e.g., Kormendy 
et al. 2009). The mass of the MBH is Mbh = 10"^M*, 
in agreement with the Magorrian et al. (1998) relation. 
The dark halo has a Navarro et al. (1997; NEW) profile 
[ph oc l/(r/r;i)(l -I- r/rh)'^, with the scale radius, and 
truncated at large radii], and a total mass Mh- For each 
(Jc, the defining parameters of the stellar mass model 
were chosen accordingly to the observational constraints 
that L B follows the Faber- Jackson relation, and that L b , 
(Tc and the effective radius Re lie on the Fundamental 
Plane of ETGs (e.g., Bernard! et al. 2003). The free pa- 
rameters defining the dark matter were chosen in agree- 
ment with the results from dynamical modeling of the 
observed motions of stars, planetary nebulae and globu- 
lar clusters at small and large radii; these indicate that 
the dark matter begins to be dynamically important at 
2~3Re (e.g., Sagha et al. 1992, Cappellari et al. 2006, 
Weijmans et al. 2009, Shen & Gebhardt 2010). This 
requires that /? = Vh/Re > 1, and TZ = Mh/M^ = 3 or 
5 (the latter value corresponding to the baryon-to-total 
mass ratio of WMAP, Komatsu et al. 2009). By solv- 
ing numerically the Jeans equations for the three mass 
components in the isotropic orbits case (e.g., Binney & 
Tremaine 1987), these choices produce M*, Mh and TZe 
(the dark-to-luminous mass ratio within Rg). Reasonable 
values of Af* /Lb = (4-10) Mq/Lb,©, and7^e = 0.2-1.0 
are obtained. The main properties of a few representa- 
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five mass models are shown in Fig. [TJ 

For a consistent comparison between observed T's and 
the characteristic temperatures derived for the mass 
models, the central stellar velocity dispersion Gc must 
be the same for observed ETGs and models. Typically, 
for nearby well observed ETGs, the value of CTc is that 
of the projected and luminosity- weighted average within 
an aperture of radius Re/8. Therefore, when defining a 
mass model, the chosen value of ac was assigned to this 
quantity. Finally, streaming motions as stellar rotation 
are not considered in these models (possible heating from 
these motions is discussed by Ciotti & Pellegrini 1996). 

5. DISCUSSION 

We investigate here the relationship between the ob- 
served T's and those expected from the various sources of 
heating (stellar motions, gravitational potential, SNIa's), 
or during the escape of the hot gas. For this purpose. 
Fig. [5] shows the run with CTc of the various temperatures 
defined in Sects. 12. 2[ 12. 3[ and [Sj together with the dis- 
tribution of the observed T values from BKF (Sect. [T]); 
for the BKF sample, the value of Uc is the luminosity 
weighted average within Re/8, and has been taken from 
SAURON studies for 12 ETGs (Kuntschner et al. 2010), 
for the remaining cases from the references in the Hyper- 
leda catalog (see Tab.[T]). 

The temperatures defined in Sects. 12.21 12. 3[ and |3] 
are mass-weighted averages, which is required when dis- 
cussing energetic aspects of the gas (e.g., the energy re- 
quired for escape as measured by < T~^^ > compared 
with the input energy from SNIa's) . When a direct com- 
parison is made with observed T's, it must be noted that 
the latter coincide with mass-weighted averages only if 
the ISM has everywhere one single temperature value; if 
the gas is multi-phase, or its temperature profile has a 
gradient, a single T value measured from the spectrum 
of the integrated emission will be close to an emission- 
weighted average (e.g., Ciotti & Pellegrini 2008, Kim 
2011). This means that, since the densest region is the 
central one, the measured T's tend to be closer to the 
central values than the mass- weighted ones. The tem- 
perature profiles observed with Chandra change contin- 
uously in shape, as the emission-weighted average T de- 
creases from <1 keV to ^ 0.3 keV: they switch from a 
flat central profile that increases outward of ~ 0.5i?e, to 
a quasi-isothermal profile, to a profile with a negative 
gradient (Diehl & Statler 2008, Nagino & Matsushita 
2009). Therefore, the lowest observed T's, presumably 
associated with the last category, may be larger than 
mass-weighted values, intermediate T's may be the clos- 
est to mass- weighted averages, while the largest observed 
T's may be lower than mass-weighted averages. An- 
other aspect to recall is that the temperatures defined 
in Sects. 12.21 12.31 and [3] refer to a gas distribution with 
pgas oc p*; this is appropriate for the continuously in- 
jected gas (e.g., for < Tinj >), while it may be less ac- 
curate when comparing observed T's with < T^^^ >, or 
when discussing the energetics of the whole gas content 
of an ETGs by means of < T^.^^ > and < T~^^ >, since 
the bulk of the hot ISM may have a different distribu- 
tion from the stars. For example, the observed X-ray 
brightness profile of gas-rich ETGs was found to follow 
the optical one, which was taken as evidence that roughly 



Pgas oc (e.g., Sarazin & White 1988, Fabbiano 1989). 
For gas-poor ETGs hosting galactic winds, the modeling 
shows that the profile pgas (f) will again be shallower than 
p*{r), though not as much as in the previous case (see, 
e.g.. White & Chevalier 1983). If pgas has a flatter radial 
profile than p^,, then it is easy to show that its mass- 
weighted < T+3^^ > will be larger than derived using 

Eq. m and its mass-weighted < T^^^^ > and < T^^J" > 
will be lower than derived using Eqs. Eland [71 In conclu- 
sion, the comparison of observed T's with mass- weighted 
expectations is the best that can be done currently, in a 
general analysis as that of the present work, though with 
the warnings above. Note, however, that all arguments 
and conclusions below remain valid or are strenghtened, 
when taking into account the above considerations about 
observed T's, or about the modifications to < Tj^^^ >, 
< Tg-,, > and < T<f,f >. 

5.1. Observed and predicted temperatures in the T — ac 

plane 

In the left panel of Fig. [51 the observed T's are com- 
pared with approximate estimates of the stellar temper- 
ature Tct, of Tinj, and of the escape temperature AT^- 
(Sect. [3]). The gas luminosity is also indicated with differ- 
ent colors, having grouped the Lx values in three ranges, 
chosen to have a roughly equal number of ETGs in each 
range. This grouping gives an indication of the gas fiow 
status, based on previous works: a galactic wind leav- 
ing the galaxy with a supersonic velocity has Lx < 10"^® 
erg s~^ (e.g., Mathews & Baker 1971, Trinchieri et al. 
2008), global subsonic outfiows and partial winds can 
reach Lx ^ 10"° erg s'^ (Ciotti et al. 1991, Pellegrini & 
Ciotti 1998), and a central infiow becomes increasingly 
more important in ETGs of increasingly larger Lx ■ Ma- 
genta ETGs (10^^ erg s^^ < Lx < 1.5 x 10^^ erg s'^) 
should then host winds, subsonic outflows, and partial 
winds with a very small inflowing region of radius < 100 
pc; cyan ETGs (1.5 x lO^^ erg s'^ < Lx < 1.2 x 10"° 
erg s~^) should host subsonic outflows and partial winds 
with an increasingly larger inflowing region (of radius up 
to a few hundreds pc); black ETGs are hot gas-rich and 
mostly inflowing 0. The best fit found for X-ray bright 
ETGs is also shown in the left panel (O 'Sullivan et al. 
2003), and gives a good representation of the distribution 
of observed T's down to a range of low temperatures and 
gas contents never explored before. The slope of the fit 
(T cx al '^^) and that of the cx cr^ relation are similar, 
with the fit being shallower; this could be due to the fact 
that not all heating sources depend on a^, see, e.g., the 
important SNIa's contribution in Tinj, that produces a 
much flatter run of T^j with tJc (Fig. [2]). The fit was 
mostly based on gas-rich ETGs, whose T's show a trend 
with CTc closer to that of T^ (an aspect further addressed 

^ In the X-ray faintest ETGs, the gas emission Lx is ~ (1 — 2) 
times the integrated emission from the population of weak unre- 
solved stellar sources (Sect.[TJ, generally referred to as AB-|-CVs. 
BKF derived an AB-I-CV emission model by jointly fitting the M31 
and M32 spectra, and then tested how the measurement of the gas 
properties may be affected by the adopted AB+CV model (their 
Sect. 3.1). By fitting with different AB+CV models, all within the 
uncertainties in the adopted one, the measured T changed only 
negligibly, for the 6 lowest-Lx ETGs of their sample; a systematic 
uncertainty of 10% - 20% was found for the gas flux (much lower 
for the X-ray brighter ETGs). 
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below in Sects. and [51) . while gas-poor ETGs depart 
most from it, since their T's change little for largely vary- 
ing (7c (as found by BKF; see Sect. 15.4] below). 

The right panel of Fig. [2] shows the stellar tempera- 
ture < >, the injection temperature < Tinj >, the 
escape temperature < T~g^^ >, and the characteristic 

temperature for slowly outflowing gas < T^^J' >, cal- 
culated for a set of representative galaxy mass models 
(Sect. [4|). At any fixed tic and Sersic index n, < T* >, 

< Tgi^av > and < T/g"'' > are larger for larger galaxy mass 
(TZ), and mass concentration (smaller The dashed 
lines represent a reasonable upper limit to the values of 
each of the characteristic temperatures, since they corre- 
spond to the most massive model ETGs, with the most 
concentrated dark matter allowed for by recent studies 
(Sect. 131). The < T~.^^ > curves lie below the simple 
approximation of the escape temperature given by -iT^; 

< T-^^ >= 4.8 < n >, for 7^ = 3, and ~ 5.2 < >, 
for the three cases with TZ = 5. 

As expected, all the < > curves lie below T^r, 
that overestimates the kinetic energy associated with the 
stellar random motions (Sect. 12.21) . Note that, from 
the virial theorem, < T* > is independent of orbital 
anisotropy, that just redistributes differently the stellar 
heating within a galaxy; the presence of ordered rotation 
in the stellar motions, instead, requires a more careful 
consideration. For any fixed galaxy mass model, this ro- 
tation would leave the total stellar heating unchanged or 
lower it, depending on whether the whole stellar stream- 
ing motion is converted into heat, or just a fraction of it 
(Ciotti & Pellegrini 1996). For the worst case that the 
stellar rotational motion is not thermalized at all, and 
the galaxy is a flat isotropic rotator, < > in Fig. [5] 
should be an overestimate of ^ 30% of the temperature 
corresponding to the stellar heating (Ciotti & Pellegrini 
1996); the possible reduction of < T* > to be considered 
should be lower than this, as far as the massive ETGs 
in Fig. [2] are less flattened and more pressure supported 
systems (e.g., Emsellem et al. 2011). 

All observed T's are located above < >; thus, addi- 
tional heating with respect to the thermalization of the 
stellar kinetic energy is needed, as noticed previously us- 
ing (e.g., Davis & White 1996, BKF). The gas could 
retain memory of its injection temperature, and have 
the additional infall heating, as examined in Sects. 15.21 
and 15.31 below. 

Finally, the values of < Tinj > ioi f — 0.85 are by far 
the largest temperatures of Fig. [21 larger than < T^,.^^ > 
up to (Tc ~ 250 km s^^; therefore, SNIa's should cause 
the escape of the gas for all ETGs up to this CTc, since the 
gas at every time is injected with an energy larger than 
required to leave the galaxy potential. This expectation 
is fulfilled by ah ETGs with CTc<200 km s"!; their X-ray 
properties (a low Lx, and T's of the order of < T^^J' >) 
agree well with what expected if outflows are important 
in them. This result had been suggested previously based 
on the low observed Lx ; now for the first time it can be 
confirmed based on the observed T values. At CTc > 200 
km s~^, instead, ETGs may have Lx far larger than 

* All the rest equal, these temperatures are also larger for smaller 
n, due to the galaxy being more massive to reproduce the same ctc, 
since the stellar mass profile is less steep (e.g., Fig.[TJ. 



expected for outflows (black symbols), and most ETGs 
where likely outflows are important (magenta or cyan 
symbols) have T much lower than < T^^J' >; these find- 
ings are discussed in Sect. 15.4] below. 

5.2. Gravitational heating in gas-rich ETGs 

We examine here the possibility that the additional 
heating with respect to the thermalization of the stellar 
kinetic energy is provided by infall heating and SNIa's. 
Davis & White (1996) assumed that in all ETGs the hot 
gas is infiowing, and suggested that the observed temper- 
atures are larger than Tcr just due to the luminous parts 
of ETGs being embedded in dark matter halos dynami- 
cally hotter than the stars; i.e., a form of "gravitational 
potential" way for the gas heating was invoked. This 
way can consist of an effect of the dark halo on the stel- 
lar motions, that are then thermalized, or directly on 
the gas during infall (e.g., via E^j.^^). The first possibil- 
ity is excluded by the < T* > curves in Fig. [21 that are 
always lower than Tg-, and that, through the Jeans equa- 
tions, include the effect of a massive dark halo consistent 
with the current knowledge of the ETGs' structure. In 
the second possibility of heating from gas infall, T^Jav is 
indeed potentially an important source of heating, that 
increases with the amount and concentration of the dark 
matter. This can be judged from Fig. [2l after consider- 
ing that < T~av >~ 5 < Th. >, and that the tempera- 
ture possibly attainable from infall was estimated to be 
T^„fl< < Tg-^v > (end of Sect. [3]). Note that Ti„fi {ii it 
behaves as < T+^y >) could he oc cr-^, a trend close to 
that shown by the T's of gas-rich ETGs (BKF; see also 
Sect. [6] below). 

Inflowing ETGs can also beneflt of the SNIa's energy 
input; for the unit mass of injected gas, this is written 
as Ey}^ = RsnEsn/M*- Both E+.^^ and E^'}^ then con- 
tribute to the required additional thermal energy with 
respect to that gained from the stellar random motions, 
i.e., to AEth = 3fc(T- < T, >)/2^mp. E+^^ and E^"^ 
can be in large part radiated in gas-rich ETGs, but they 
seem to far exceed the required AEth- For example, 
for the highest-Lx of Fig. H AEth (1 - 2) x lO*® 
erg Mq^ (i.e., ~0.2-0.4 keV), when adopting an aver- 
age galaxy mass model as that of the thick black line in 
Fig. [21 The energy spent in radiation can be estimated, 
in a stationary situation, as Lx /M* (per unit injected 

gas mass); using Lx from BKF and deriving Af* as in 
Sect. 12. ll for the same distances in BKF and galactic B- 
magnitudcs given by Hyperleda, for the gas-rich ETGs 
one flnds that Lx/ M^, ranges between (0.5 — 3.4) x 10''® 
erg Mq^ . The energy available is far larger than the 
sum of AEth and Lx/M*: E]?^ = 7.3 x 1048(//0.85) 
erg Mq^, and 0.5i?+av ranges from 4 x 10*® erg Mq^ 
{(Jc ^ 220 km s-i) to 8 X erg Mq^ (ctc - 300 km 
s~^), for the mass model with the thick black line in 
Fig. [21 These results are detailed in Fig. [3l where the 
values of AEth and Lx/M^, for each galaxy are shown, 
together with various combinations of i^^^v ^^'^ ^s'n- 
conclusion, additional input energy for the gas to account 
for the observed T's of gas-rich ETGs seems available in 
a sufficient amount, even if / were to be < 0.85. 
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5.3. Outflows and SNIa's heating 

The gas is not mostly inflowing in all ETGs, while it 
is hotter than < T* > in all of them. When in outflow, 
the radiative losses are far smaller, but energy is spent in 
extracting the gas from the galaxy and giving it a bulk 
velocity. We discuss here the possibility of heating from 
the SNIa's energy input to account for the observed AEth 
of ETGs where outflows are likely important (those with 
low/medium Lx, magenta and cyan symbols in Fig. [2|). 

We assume that the SNIa's energy is used for the uphft 
of the gas and the kinetic energy with which it escapes 
from the galaxy, and, neglecting radiative losses, that all 
the remaining part is available to account for the ob- 
served AEth- Then, the energy balance per unit mass 
of injected gas is E^'^ = E'^^ + Eout + AEth, where 
AEth is the same as in Sect. 15.21 and Eout — is 
the mass-averaged kinetic energy of the escaping mate- 
rial per unit gas mass. Figure |4] shows AEth derived 
from this balance, for Vout = Cs, and calculated for 
7 = 5/3 and kT = 0.3 keV, a temperature of the or- 
der of that observed for ETGs likely in outflow (Fig. [2]). 
Adopting Vout ~ Cs (independent of ctc) produces an Eout 
on the upper end of those expectetfl, and then the esti- 
mate of AEth may be biased low. With Vout ~ Cg, Eout 
is just ~ 0.1(//0.85)£'c.°^. The energy needed for gas 
extraction E~j.^^, for the same galaxy mass model used 
in Sect. [El varies instead from ~ l/3(//0.85)£;|°^ for 
ac = 150 km s~\ to - Ef^ for CTc = 250 km s'^; this 
explains the strong dependence of the predicted AEth on 
(Tc in Fig. m It is clear from this figure that for / = 0.85 
SNIa's can account for the needed heating in all ETGs 
with low/medium Lx] for / = 0.35, instead, the temper- 
ature increase would fall short of what required for all 
ETGs. In Fig. HI the energy losses due to radiation are 
also shown; their small size supports the hypothesis that 
in most cases they do not affect significantly the energy 
budget of the gas. 

Given the flat distribution of the observed points in 
Fig. Ill and the steep behavior of the curves predicting 
AEth, the value of / required to account for the ob- 
served AEth increases with Cc- In particular, the value 
of / ^ 0.85 that is required at high CTc would produce 
an expected AEth at low cFc that is larger than ob- 
served. A possible solution could reside in the efficiency 
of the SNIa's energy mixing process. In massive, gas-rich 
ETGs, SNIa's bubbles should disrupt and share their en- 
ergy with the local gas within ^ 3 x 10^ yr (Mathews 
1990); for a Milky Way-size bulge in a global wind, in- 
stead, 3D hydrodynamical simulations of discrete heating 
from SNIa's suggest a non-uniform thermalization of the 
SNIa's energy, with overheated gas by a SNIa explosion 
at the bulge center that is advected outwards, carrying 
a large fraction of the SNIa energy with it (Tang et al. 
2009). For subsonic outflows the mixing is expected to 
be more local and more complete (Lu & Wang 2011). 
The magenta and cyan ETGs in Fig. |4] have gas den- 
sities and luminosities larger than those considered by 

^ For example, in a wind solution, the terminal (i.e., the largest) 
velocity of the flow is roughly the central sound speed (White & 
Chevalier 1983); moreover, in this solution the gas is likely to be 
already too "fast" with respect to that of most magenta ETGs in 
Fig. El due to their Lx (e.g., Trinchieri et al. 2008). 



Tang et al. (2009); however, if a discrete heating effect 
were still present at CTc < 200 km s~^, it could quali- 
tatively explain a lower / for these galaxies. There is 
also the possibility that < > has been overestimated 
(and then the observed AEth underestimated) at the low 
(Tc, if these ETGs are less pressure supported systems 
(e.g., Emsellem et al. 2011), and the stellar rotational 
streaming is not all thermalized (Sect. 15. ip . Another 
possible explanation could be that ETGs with CTc > 200 
km s^^ are less outflow-dominated than those at lower CTc 
(though this is not supported just based on Lx , since ma- 
genta ETGs are found over the whole CTc range in Fig.|4]), 
so that their iiwav would be lower than assumed by the 
curves in Fig. |4l and more SNIa's energy would be avail- 
able for heating. In this way, / could have a value < 0.85, 
possibly similar for all ETGs. 

Finally, a comparison of Figs. [3] and |3] shows that the 
average AEth is slightly larger for the X-ray brightest 
ETGs (for which it ranges between 0.1-0.5 keV) than for 
the X-ray faintest ones (0-0.3 keV); moreover, while the 
AEth in Fig.[3]can be explained even with / < 0.85, it is 
required that / ~ 0.85 for the X-ray faintest ETGs with 
the largest ac in Fig. |31 Both facts are the consequence 
of the large fraction of the SNIa's energy input that is 
used in gas extraction where outflows dominate, while 
all the SNIa's energy remains within the galaxies where 
inflow dominates. 

In conclusion, even for ETGs with low/medium Lx, 
a fundamental X-ray property as T can be accounted 
for by simple arguments, just based on realistic galaxy 
mass models, and reasonable SNIa's heating capabilities. 
There may be, though, more energy available for the gas 
in ETGs with Uc < 200 km s~^ than can be accounted 
for by the present simple scenario. 

5.4. The temperature and gas flows status in ETGs of 
intermediate mass 

The observed X-ray properties (low Lx and T ^< 
'^csc >)' and the energy budget of the gas (e.g., < 
^gmv > VS. < T,nj >), for ETGs with CTc<200 km s'^ 
are all consistent with the expectations for outflows; the 
large Lx and < T^.^^ > larger than < Tinj > of ETGs 
with (Tc > 250 km s^^ agree with the gas being mostly 
inflowing. For 200 < (Tc(km s^^) < 250, instead, ETGs 
show very different Lx and T, whose values seem unre- 
lated to < T^^^ > and to the relative size of < T~^^ > 
and < Tinj > (Fig. [2]). For example, for / = 0.85, 
< Tinj > exceeds < T~^^ >, but most T values lie well 

below < r^'^^ >, and even high Lx values (incompatible 
with outflows) are common. One first explanation could 
be that / < 0.85; for example, for / = 0.35, < T^j > be- 
comes lower than < T~^^ > at ac ^ 180 km s^^ (Fig. [2]). 
Four ETGs with ac > 200 km s~^ and a very low Lx 
(Fig.d magenta symbols), though, require that / > 0.35 
in them, and then that / varies from galaxy to galaxy, or 
that their gas was removed by other processes as an AGN 
outburst (e.g., Machacek et al. 2006, Ciotti et al. 2010), 
or a merging or an interaction (Read & Ponman 1998, 
Sansom et al. 2006, Brassington et al. 2007). An event 
like the latter two in the recent past is unlikely for three 
of these ETGs (NGC1023, NGC3115, NGC3379), that 
are very regular in their stellar morphological and kine- 
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matic properties , while in the other ETG (NGC4621) 
that hosts a counter-rotating core (Wernli et al. 2002) it 
is possible. 

A second explanation could be that < Tinj > exceeds 

< Tgi^av > only at the present epoch: while < T, > and 

< T~^^ > are independent of time, Tsn may have been 
lower in the past (Eq. [T]), to the point that < Tinj > 
may have been lower than < T~.^^ > for more ETGs 
than in Fig. [2] (that represents a snapshot of the present 
epoch). The gas then could have accumulated and ra- 
diative losses have become important, even for the gas 
injected in later epochs. In fact, the population synthe- 
sis models of Sect. [2] predict that was larger at early 
times (e.g., by ~ 6 times at an age of 3 Gyr), and then 
to keep Tsn high in the past, from Eq. [1] Msn must 
decrease with time i at a rate similar to or steeper than 
that of the stellar mass losses (M, cx t~^-^; Ciotti et al. 
1991). Recent observational estimates indicate instead 
a SNIa's rate decaying close to (Maoz et al. 2010, 
Sharon et al. 2010), thus Tsn and < T^j > should be 
increasing with time, reaching the values of Fig. [2] at the 
present epoch. Then, a "cooling effect of the past" would 
explain a moderate or high Lx even where < Tnj > ex- 
ceeds < T-^^ > in Fig. O 

Both explanations, / < 0.85 and/or a lower < Tmj > 
in the past, can account for the lack of a widespread 
presence of outflows for 200 < crc(km s'^) < 250. Both, 
though, require a mechanism different from the SNIa's 
energy input to cause degassing in some low Lx ETGs 
in the same Cc range. 

Another possibility is that partial winds become com- 
mon for cTc > 200 km s~^: these ETGs host an inner 
inflow and an outer outflow (e.g., MacDonald & Bailey 
1978), with variations in the galactic structure causing 
different sizes for the inflowing region, and then differ- 
ent Lx (Pellegrini & Ciotti 1998). This possibility holds 
even for / = 0.85, and for both kinds of time evolu- 
tion of < Tinj >; in fact, if the flow is decoupled, ETGs 
may host a central inflow even if < Tinj > is larger than 

< Tg^av > ■ Similarly, the observed T's can be lower than 

< T^^ >, in ETGs where the outflow is only external; 
in this case, the observed T's may also be lowered by 
radiative losses in the central inflowing region. 

6. THE TEMPERATURE OF INFLOWS AND OUTFLOWS 

In Fig. [3 the hottest gas is in ETGs with the highest 
Lx, and the coolest one in ETGs with the lowest Lx 
(as also found by BKF). This feature is also present in 
the Lx - (Jc relation, where ETGs with kT > 0.4 keV 
are the X-ray brightest (with one exception), while those 
with kT < 0.3 keV are the X-ray faintest (BKF). AU 
this may seem contrary to the simple expectation that 
hotter gas is needed for escape, and that the hotter the 
gas, the stronger the outflow, the lower the gas content. 
We re-examine below this point, first across the whole (Tc 
range, and then at fixed CTc- 

A proper consideration of whether outflowing ETGs 
possess hotter or colder gas than inflowing ones requires 
that all T's are rescaled by a temperature equivalent to 
the depth of the potential where the gas resides (for ex- 
ample, by Tcr). Is there a trend then of the distance of 
the observed T's from < T/.f >, or < > ? This is 



examined by Fig. [Sj where temperatures are rescaled by 
Tcr, and the dependency of all the curves in Fig. [2] 
is removed. Figure [S] shows that for CTc < 200 km s~^ 
the observed points reach < T^^ > /To-, a result similar 
to that of Fig. [2l and that they fall below it with in- 
creasing (Tc (with a transition region of large dispersion 
in T/T^). Therefore ETGs with ctc < 200 km s"! are 
indeed the hottest, relatively to the virial temperature; 
since in these ETGs outflows are important (Sect. lOT) . 
indeed the flow is relatively hotter in outflows, and T in- 
creases going from outflows to inflows only in an absolute 
sense. The T's of the X-ray brightest ETGs should show 
a dependence on a^, if gravitational heating of the gas 
dominates over SNIa's heating, and then they should lie 
within a horizontal zone in Fig. [51 The observed distri- 
bution does not disagree with this kind of dependence, 
but more cases are needed to flrmly establish its pres- 
ence; such a dependence is not expected, though, since 
the SNIa's heating should easily dominate over the grav- 
itational one (Sect. 15. 2p . 

We finally compare the T's at similar tJc, in the most 
populated region of Fig. [SJ for 200 < a-c(km s^^) < 250. 
Here the variation of T/Tg. is the largest, and is covered 
by ETGs of all X-ray emission levels; the X-ray bright- 
est ETGs are found at T/T^ > 1.1, while the lowest 
T/Tcr values belong to the X-ray fainter ETGs. For ex- 
ample, two of the three lowest T/To- values of the fig- 
ure (those of NGC3379 and NGC4621) belong to the 
X-ray faintest group. While heating sources seem abun- 
dant in gas-rich ETGs to account for their T's (see, e.g.. 
Fig. [31 and the additional possibility of MBH heating. 
Sect. 12.41) . even after taking into account their radiative 
losses, this result remains more difficult to explain for 
ETGs of low/medium Lx, and may require ad hoc so- 
lutions. It may be another representation of what men- 
tioned in Sect. 15. 3[ that AEth can be larger for the X-ray 
brightest ETGs than for the X-ray faintest ones, due to 
the different employment of the SNIa's input energy; or 
it could be that /< 0.85 in these ETGs so that SNIa's 
cannot make their gas hotter than this (Fig. [J) ; or their 
galaxy structure may be much different from an average 
one, so that dividing all T's for the same To- produces a 
biased view; or the evolutionary history of the gas may 
have been peculiar. Certainly, this trend needs further 
investigation and, if confirmed, it will provide the basis 
for further theoretical work. 

7. CONCLUSIONS 

This work has focussed on the origin of the hot gas 
temperatures recently derived for a sample of ETGs ob- 
served with Chandra down to galaxy masses and X-ray 
luminosities smaller than ever before. A few character- 
istic mass-weighted average temperatures have been de- 
fined for a gas distribution pgas (r) oc p, (r) , as for the gas 
shed by stars: the virial temperature < T, >; the injec- 
tion temperature < Tinj >, as the sum of < T, > and 
of a temperature equivalent to the SNIa's kinetic energy 
input (with a factor / allowing for its uncertain thermal- 
ization); the escape temperature < T~av >, defined as 
the temperature equivalent of the energy required for es- 
cape from the gravitational potential; a fiducial value for 
the temperature of escaping gas, evaluated on a stream- 
line of very subsonic velocity (< T^^J' >= 0.6 < Tg^^v >); 
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and finally, the temperature equivalent to the energy lib- 
erated by the gas inflow to the galactic center, < T^^^ >. 
These temperatures were then calculated for a set of rep- 
resentative galaxy mass models, made by the superposi- 
tion of a central MBH, and a stellar and a dark mass 
density distributions, with parameters constrained from 
the fundamental scaling laws of ETGs and recent obser- 
vational findings. The main properties of the character- 
istic temperatures are that: 

• All temperatures scale as (except for < Tinj >), 
and increase for larger and/or more concentrated mass 
content. For the adopted set of representative galaxy 
mass models, < T, > is lower than To- (by ^ 0.1 — 0.2 
keV), < Tg-^, >« 5 < T, >, and < T+^, >« 2 < 
^grav tti^ temperature that can be produced by infall 
heating, though, will be much lower than < T~^^ >, 
due to energy losses in radiation, kinetic energy of mass 
condensations, and mass drop-outs from the fiow. 

• < Tinj > is by far the largest of the characteristic 
temperatures, due to the important SNIa's contribution 
(independent of <7c); for / = 0.85, it is larger than the 
minimum injection temperature for global escape up to 
(Tc ~ 250 km s-^ 

The comparison of the characteristic temperatures 
with those observed, in the T — Uc plane, shows that: 

• The best fit T — Cc relation previously found for X- 
ray bright ETGs reproduces the average trend of the ob- 
served T down to low temperatures, and low Lx- ETGs 
with low/medium Lx show the largest departures from 
this fit, which can be explained by the variety of gas fiow 
phases possible in them (winds, subsonic outflows, par- 
tial winds), where the main input energies (from SNIa's 
and gas infall) are used in different ways. 

• All observed T's are larger than < T, >; the ad- 
ditional heating of the gas AEth, with respect to that 
provided by the thermalization of the stellar motions, is 
AEth « O'- 0.3 kcV for the X-ray faintest ETGs, and 
AEtfi « 0.1 — 0.5 keV for the X-ray brightest (for a rep- 
resentative galaxy mass model). 

• In a stationary situation, AEth of the X-ray bright- 
est ETGs can be accounted for by the energy input of 
SNIa's and gas infall, even if they are much reduced with 
respect to standard assumptions (i.e., / can be < 0.85). 
The gravitational heating produces a T oc trend, that 
may be present in the X-ray brightest ETGs; the SNIa's 
heating, though, is expected to be dominant. 

• AEth can be provided by SNIa's in X-ray fainter 
ETGs, where outflows are important; most of the SNIa's 
energy is needed for gas extraction, and less for the ki- 
netic energy of the escape. The value of / to account 
for the observed AEth increases with CTc, until the whole 
SNIa's energy (/ « 0.85) is required at the highest Cc- 
With this /, though, at low CTc the observed AEth are 
lower than expected. Possible solutions require a dif- 
ferent efficiency of the SNIa's energy mixing process, or 
an overestimate of < T, > at low Uc if these ETGs are 
less pressure supported systems, or a more complex fiow 
status than in the simple scheme adopted. 

• At low (TcS5200 km s~^, < Tinj > is larger than 
< Tg^av >j tiic Tx values are low and the T's are 
of the order of < T^i^ >: all this agrees well with 
what expected for outfiows. At high CTc > 250 km s^^. 



< Tinj > is lower than < T^^^v ^^-^ the high Lx and 
T can be explained with the gas mostly inflowing. For 
200 < tTc(km s"-'^) < 250, instead, there is a large varia- 
tion in Lx and T. Possible explanations could be that 
the SNIa's energy input varies from galaxy to galaxy, 
and/or that < Tinj > was lower in the past, due to the 
different time evolution of the mass loss and the SNIa's 
rate; or that partial winds become common, with the flow 
status less related to the values of < Tinj >, < T'g^av >i 
and < T^^J> >. 

• When measured relatively to the depth of the poten- 
tial well, the observed temperatures T/T^ are larger for 
cTc < 200 km s~^ (outflows), and lower for CTc > 250 km 
s^^ (inflows). The observed T's then increase from out- 
flows to inflows only in an absolute sense, and the gas is 
relatively hotter in outflows. In the intermediate region 
of 200 < (Tc (km s^^) < 250, lower Lx values tend to cor- 
respond to lower T and T/Tg-, which requires ad hoc ex- 
planations, and then deserves further observational and 
theoretical investigation. 

I thank Luca Ciotti for helpful discussions, and Dong- 
Woo Kim and the referee for useful comments. 
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Table 1 

Observed properties of the ETG sample. 



Name 


kT 






error 




Lx 




CTc 


Ref 




(keV) 










(10"" erg s- 




(km s^l) 




(1) 


(2) 




(3) 






(4) 




(5) 


(6) 


NGC 720 


0.54 


-0 


01 


+0 


01 




5.06 




241 


Binney et al. 1990 


NGC 821 


0.15 


-0 


05 


+0 


85 


2 


13 X 10" 


3 


200 


Kuiitschiicr ct a,l. 2010 


NGC1023 


0.32 


-0 


01 


+0 


02 


g 


25 X 10" 


2 


204 


Kuiitsctiricr ct 9.1. 2010 


i V7 V_' _L W (J 


0.34 


-0 


02 


+0 


02 


4 


37 X 10~ 


1 


215 




NGC1316 


0.60 


-0 


01 


+0 


01 




5.35 




230 


D'Onofrio ct al. 1995 


NGC1427 


0.38 


-0 


11 


+0 


26 


5 


94 X 10" 


2 


171 


D'Onofrio ct al. 1995 


NGC154q 


0.35 


-0 


04 


+0 


04 


3 


08 X 10" 


1 


210 


Loiigo ct al. 1994 


NGC2434 


0.52 


-0 


05 


+0 


04 


y 


56 X 10" 


1 


205 


Loiigo ct al. 1994 


NGC2768 


0.34 


-0 


01 


+0 


01 




1.26 




205 


Kiint^rhnpr pt al 2010 


NGC3115 


0.44 


-0 


10 


+0 


16 


2 


51 X 10~ 


-2 


239 


Fisher 1997 


NGG3377 


0.22 


-0 


07 


+0 


12 


\ 


17 X 10~ 


-2 


144 




NGG337Q 


0.25 


-0 


02 


+0 


03 


4 


69 X 10~ 


-2 


216 


Knntsirhnpr pf al 9010 


NGC3384 


0.25 


-0 


15 


+0 


17 


3 


50 X 10~ 


-2 


161 


Kimtsrhner et al 2010 


NGC3585 


0.36 


-0 


05 


+0 


06 


\ 


47 X 10~ 


1 


198 


Fisher 1997 


NGC3Q23 


0.45 


-0 


01 


+0 


01 




4.41 




250 


Pellea-rini et al 1997 


NGC4125 


0.41 


-0 


01 


+0 


01 




3.18 




227 


Bender et al. 1994 


NGr;42fi1 


0.66 


-0 


01 


+0 


01 




7.02 




300 


RpTirlor of al 1 QQ4. 


NGG4278 


0.32 


-0 


01 


+0 


01 


2 


63 X 10" 


1 


252 




NGC43fi5 


0.44 


-0 


02 


+0 


02 


5 


12 X 10" 


1 


245 


Bender et al. 1994 


NGC4374 


0.63 


-0 


01 


+0 


01 




5.95 




292 


Kuntschncr et al. 2010 


NGC4382 


0.40 


-0 


01 


+0 


01 




1.19 




187 


Kuntschncr et al. 2010 


NGC4472 


0.80 


-0 


00 


+0 


00 




18.9 




294 


Bender et al. 1994 


NGC4473 


0.35 


-0 


03 


+0 


05 


1 


85 X 10" 


1 


192 


Kuntschner et al. 2010 


NGC4526 


0.33 


-0 


01 


+0 


02 


3 


28 X 10- 


1 


232 


Kuntschner et al. 2010 


NGC4552 


0.52 


-0 


01 


+0 


01 




2.31 




268 


Kuntschner et al. 2010 


NGC4621 


0.27 


-0 


09 


+0 


13 


6 


08 X lo- 


2 


225 


Kuntschner et al. 2010 


NGC4649 


0.77 


-0 


00 


+0 


00 




ll. 7 




315 


Bender et al. 1994 


NGC4697 


0.33 


-0 


01 


+0 


01 


1 


91 X 10- 


1 


174 


Binney et al. 1990 


NGC5866 


0.35 


-0 


02 


+0 


03 


2 


42 X 10- 


-1 


159 


Fisher 1997 



Column (1): galaxy name. Cols. (2), (3) and (4): the hot gas temperature, its uncertainty, and the 0.3—8 keV gas luminosity, from BKF. 
Col. (5): the stellar velocity dispersion, as the luminosity-weighted average within an aperture of radius Re/8, with its reference in Col. 
(6). 
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R(pc) 



Figure 1. Mass (up), B-band surface brightness (middle), and projected velocity dispersion (bottom) profiles of three-component galaxy 
models (MBH+stars+dark matter), for two representative ETGs with isotropic orbits and an aperture velocity dispersion within Re/ 8 of 
Uc = 260 km s^^ (solid lines; Lb = 5 x IO^'^Ls^q, and Re = 6.5 kpc), and of ac = 200 km s""*^ (dashed lines; Lb = 2 x 1O^'^-Ls^0, and 
Re = 3.6 kpc). Red lines refer to a stellar Sersic profile with index n = 4, blue ones with n = 5. The dark halo in the upper panel (black, 
with the same line type as the corresponding stellar profile) follows the NFW profile, with /3 = 2, and TZ = 3 (for n = 4) or 7^ = 5 (for 
n = 5), and then 7?.e=0.24 or 0.41, from the Jeans equations (see Sect. |4]l. 
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a_, (km s-') CT_, (km s^i) 



Figure 2. The relationship between the observed gas temperature (from BKF) and ac (see Sect.[5TT}. Symbols surrounded by a magenta 
and cyan circle have respectively 10"^* erg < Lx < 1-5 X 10"^^ erg s~^, and 1.5 X 10^^ erg < Lx < 1-2 X 10'**' erg s~^; all other 
ETGs have larger Lx- Left panel: in green To- (Sect.O, and the simple estimate of 4To- for the escape temperature (Sect.[3ll; in blue the 
best fit o-c oc T°-56±0 09 found from ROSAT data (O'SuUivan et al. 2003); in black two cases of T^^j (Eq. El, calculated using T^. Right 
panel: 1) < T, > (Eq.[2l lowest bundle of lines), calculated for four representative galaxy mass models (made of MBH+stars+dark halo), 
with a Sersic index n = 4 (black lines) or n = 5 (red line), and the dark matter parameters TZ., f), Tie indicated on each curve (Sect. [Jjl; 2) 
< r^av > and < Tgj,"'' > (Sect.[3| for the same mass models adopted for < T, >, with the corresponding line type and color; 3) < Ti„j > 
calculated using < T« >, with the corresponding line type and color, and / = 0.85 or / = 0.35 in Eq. |3] 
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150 200 250 300 

CT^ (km s"') 



Figure 3. The run with Cc of the energies provided by SNIa's {Ej?^) and gas infall (-B^av); points with errorbars show, for the ETGs 
with the largest Lx in Fig. [2] the additional thermal energy {AEff^) with respect to that gained from the thcrmalization of the stellar 
random motions, required to explain the observed T's (see Sect. [5T2t ; all these energies {E^'^, B^avi a^nd AE^f^), that are defined per 
unit mass in the text, have been multiplied by 2fimp/3 to obtain their temperature-equivalent in keV plotted here. Solid lines show the 
sum of the energies provided by the SNIa's and infall, for two cases of /, and for E^av rescaled by a factor of 0.1, and calculated for the 
galaxy mass model described by the thick black line in Fig. [2] (right panel). The dashed line gives for reference the value of 0.1-E^avi and 
the dot-dashed line the value of Eg'^ for / = 0.35. The points with errorbars are obtained subtracting to the observed kT the value of 
A: < T* > corresponding to its cTc, for the mass model adopted for -E^av; each point is linked by a red line to a point including the energy 
spent in radiation observed for that ETG (i.e., the upper point measures AEtf^ + Lx/M^). See Sect. IS^ for more details. 
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Figure 4. The run with of the energy provided by SNIa's (E^g^, for three cases of /), after subtraction of the energy needed for the 
removal of the gas from the galaxy (-E^av)i and for escape (Eout) at an average Vout = Cs(0.3 keV); all energies have been computed as in 
Fig. [3] The adopted galaxy mass model is that corresponding to the thick black line in Fig. [2] (right panel). Points with errorbars show the 
additional thermal energy required to explain the observed T's, calculated as in Fig. [S] for ETGs with low/medium Lx in Fig. [2] For the 
cyan ETGs, a red line connects each point with the value including the radiated energy, as in FigJS] the red lines of the magenta ETGs, 
whose Lx values are the lowest, would be included within the colored circle, if shown. See Sect. l5.3l for more details. 
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150 200 250 300 



(km s"') 

Figure 5. The same as in Fig. [2] (riglit panel), witli temperature values rescaled by To-. 



